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ABSTRACT 

The statistical validation of transiting exoplanets proved to be an efficient technique to secure 
the nature of small exoplanet signals which cannot be established by purely spectroscopic 
means. However, the spectroscopic diagnoses are providing us with useful constraints on the 
presence of blended stellar contaminants. In this paper, we present how a contaminating star 
affects the measurements of the various spectroscopic diagnoses as function of the param¬ 
eters of the target and contaminating stars using the model implemented into the FASTIS 
planet-validation software. We find particular cases for which a blend might produce a large 
radial velocity signal but no bisector variation. It might also produce a bisector variation anti¬ 
correlated with the radial velocity one, as in the case of stellar spots. In those cases, the full 
width half maximum variation provides complementary constraints. These results can be used 
to constrain blend scenarios for transiting planet candidates or radial velocity planets. We re¬ 
view all the spectroscopic diagnoses reported in the literature so far, especially the ones to 
monitor the line asymmetry. We estimate their uncertainty and compare their sensitivity to 
blends. Based on that, we recommend the use of BiGauss which is the most sensitive diagno¬ 
sis to monitor line-profile asymmetry. In this paper, we also investigate the sensitivity of the 
radial velocities to constrain blend scenarios and develop a formalism to estimate the level of 
dilution of a blended signal. Finally, we apply our blend model to re-analyse the spectroscopic 
diagnoses of HD16702, an unresolved face-on binary which exhibits bisector variations. 

Key words: techniques: spectroscopic ; techniques: radial velocities ; methods: data analysis 
; planetary systems ; binaries: spectroscopic 


1 INTRODUCTION 

One of the main objectives of detecting new extrasolar planets is 
to constrain theories of planet formation, migration and evolution 
(e.g. Mordasini et al. 2009). Exoplanet discoveries are therefore 
used to build up our general knowledge about planetary systems. 
It is therefore crucial that these exoplanets are indeed bona-fide 
planets. However, periodic signals with a relatively low amplitude 
compatible with a planetary nature are not necessarily planets 
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(Cameron 2012). Indeed, many astrophysical non-planetary 
scenarios (the so-called “false-positives”) can mimic the signal of 
a planet. All techniques are affected: the radial-velocity planetary 
signals might be mimicked by stellar activity (e.g. Queloz et al. 
2001) or nearly face-on binaries (Diaz et al. 2012; Wright et al. 
2013) ; the photometric transit signals of a planet can be mimicked 
by various configurations of eclipsing binaries (e.g. Morton & 
Johnson 2011; Santeme et al. 2013a) ; false-positive scenarios 
might also produce misinterpreted direct imaging and microlens- 
ing exoplanet detections (Zurlo et al. 2013; Hwang et al. 2013, 
respectively). In some cases, the source of the signal is a planet, but 
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the amplitude of the signal is diluted by the presence of a brighter 
star in the system (e.g. Santos et al. 2002). Therefore, the derived 
properties of the planets (without accounting for the contaminating 
light) might be significantly different, hence limiting the accuracy 
of statistical analysis based on planet detections (e.g Howard et al. 
2012; Petigura et al. 2013; Fressin et al. 2013). 

It is therefore fundamental to secure the planetary nature of 
a detection. This is especially important for transiting exoplanets, 
for which a substantial fraction of transit-like events are not 
of planetary origin (e.g. Almenara et al. 2009; Moutou et al. 
2009; Santerne et al. 2012b; Coughlin et al. 2014). One way to 
establish the planetary nature of a transiting candidate is through 
high-resolution spectroscopy and radial velocity measurements 
(Konacki et al. 2003). Spectroscopic observations can efficiently 
screen out eclipsing binaries as well as triple systems (Santerne et 
al. 2012b). Then, to establish the planetary nature of a candidate, 
it is commonly considered that a transit candidate has to satisfy 
all the following requirements: (1) a radial velocity variation is 
significantly detected in phase with the transit ephemeris ; (2) the 
inferred mass is compatible with that of a planet; (3) no significant 
correlation is found between this RV variation and the asymmetry 
of the line profile (the bisector) ; (4) the signal is stable in time and 
amplitude (Santos et al. 2014). 

This technique works relatively well for systems where the 
reflex motion of the host can be observed with current instrumen¬ 
tation. However, for the thousands of potential small planets that 
the Kepler space telescope detected (e.g. Batalha et al. 2013), the 
reflex motion is too small to be observed (Santerne et al. 2013b), 
therefore limiting the establishment of their planetary nature. 
Other non-spectroscopic methods exist to measure the mass of a 
transiting planet, such as the relativistic Doppler-boosting effect 
(Mazeh & Faigler 2010) and through dynamical interaction with 
other planets in the system (Holman et al. 2010). However, these 
techniques are limited to very short orbital period planets or 
multiple-planet systems (Shporer et al. 2011; Boue et al. 2012, 
respectively). 

An alternative solution has been proposed by Torres et al. 
(2011), called the planet-validation. This technique consists in 
evaluating the probability of each scenario to produce the observed 
data. If the planet scenario turns out to be significantly the most 
likely scenario, then the planet is considered as validated statisti¬ 
cally. However, as demonstrated by Diaz et al. (2014), if the transit 
signal has a signal-to-noise ratio below 150, it is not possible to 
validate the nature of a transiting exoplanet based only on the 
photometric data. In such configuration, the validation relies on the 
accuracy of the a priori information. This is clearly a limitation 
for the exploration of planets like the Earth since the statistics in 
this population is extremely low. To overcome this limitation, the 
solution is to increase the statistical weight of the data compared 
with the priors in the validation procedure, by modelling data from 
different techniques providing different constraints to the various 
scenarios. 

The objective of this paper is to show how spectroscopic 
diagnoses can be used to constraints false-positive scenarios in the 
context of the validation of planets performed with the FASTIS 


software* (Diaz et al. 2014). For that, we explore the effects of 
a blended star on the measurement of both the radial velocity, 
the full width half maximum (FWHM) and the line-asymmetry 
diagnoses already defined in the literature (e.g. the Bisector 
Inverse Slope). These diagnoses are commonly used to disentangle 
planetary signals from stellar activity in radial velocity data (see 
e.g. Dumusque et al. 2012, and reference therein). These diagnoses 
have already been used by, e.g. Santos et al. (2002), Torres et 
al. (2004), Torres et al. (2005), Dfaz et al. (2012), Wright et al. 
(2013), and recently by Santerne et al. (2014) to analyse blended 
radial velocity of double-line binary systems mimicking planets or 
brown dwarfs. 

Note that the impact of stellar contamination on precise radial 
velocity has already been studied by Cunha et al. (2013) in the con¬ 
text of the preparation of the ESPRESSO spectrograph. Previously, 
the effect of a contaminating star has also been studied by Torres et 
al. (2005) to constrain the false-positive scenario of OGLE-TR-56 
which mimic both a transit and radial velocity signals of a planet. 
In this paper, we extend the work done by these authors, including 
all the spectroscopic diagnoses. We also investigate how these 
diagnoses change as function of the parameters of the target and 
contaminating star. Einally, we propose a formalism to estimate the 
amplitude of a radial velocity signal as function of the parameters 
of the blended components. 

We first explore in Section 2 the impact of a blended con¬ 
taminant on precise radial velocity and the line-profile diagnoses. 
Then, we present in Section 3 the sensitivity of radial velocity data 
to constrain blend scenarios. In Section 4, we test our FASTIS 
RV-blend model by re-analysing an unresolved nearly face-on 
binary, HD16702. We report some caveats in Section 5. Finally, 
we draw our conclusions in Section 6 and discuss how these 
radial velocity constraints complement the photometric ones in the 
context of planet validation. 

Several appendices are available at the end of this paper. In 
Appendix A, we review how to compute the various spectroscopic 
diagnoses used in this paper. In Appendix B, we describe our 
method to simulate the line profile of a blended stellar system. In 
Appendix C, we present the computation time as well as the numer¬ 
ical precision and accuracy to measure each line-profile diagnosis. 
Finally, in Appendix D, we provide the parameters used to analyse 
the system HD16702 with FASTIS, together with their prior and 
posterior distributions. 


2 EFFECT OF A CONTAMINATING STAR ON THE 
SPECTROSCOPIC DIAGNOSES 

2.1 Presentation of the spectroscopic diagnoses 

At the resolution of most spectrographs used to measure precise 
radial velocities of stars (40 000 - 120 000), an unblended and un¬ 
saturated weak line profile of a slow-rotating star could be approx¬ 
imated by a Gaussian function. Averaging thousand of such lines 
using the cross-correlation technique (hereafter CCF Baranne et al. 
1996; Pepe et al. 2002) also produces a Gaussian profile. The CCF 
presents the advantage of having a much higher signal-to-noise ra¬ 
tio compared with individual lines without modifying the shape of 
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the averaged line profile. This function is characterised by a given 
location in the radial velocity space, a contrast and a FWHM. An 
asymmetry in the observed line profile can be produced by, e.g. 
the convective blue shift, pulsations, stellar spots (Dumusque et 
al. 2014), and contaminating stars. To monitor the line asymme¬ 
try, several diagnoses have been defined in the literature. They are 
detailed in the Appendix A of this paper and illustrated in Fig. Al. 
For clarity, we list them here: 

• BIS (Queloz et al. 2001): it consists in computing the velocity 
span between the average of the top portion of the bisector and the 
bottom portion of the bisector, within some defined limits. 

• Vspan (Boisse et al. 2011a): it consists in fitting a Gaussian 
profile to the top and bottom part of the line profile. The asymmetry 
diagnosis is defined as the difference between the two velocities. 

• BiGauss (Nardetto et al. 2006): it consists in fitting an asym¬ 
metric Gaussian profile to the line profile (see eq. A7). The asym¬ 
metry diagnosis is defined as the difference between the radial ve¬ 
locity found with the symmetric and asymmetric Gaussian profiles. 

• Vasy (Figueira et al. 2013): it consists in computing the dif¬ 
ference of spectral information (as defined by Bouchy et al. 2001) 
contained in the blue and red wings of the line profile. 

• Wspan (this work, see Appendix A): it consists in fitting a 
Gaussian profile to the red and blue wings of the line profile. The 
asymmetry diagnosis is defined as the difference between the two 
velocities. 

These diagnoses can be measured directly on each individ¬ 
ual stellar lines available in the spectrum, as long as their profiles 
have been correctly reproduced by the spectrograph and are not per¬ 
turbed by an external source like other unresolved lines, or the spec¬ 
trum of a calibration lamp, an iodine cell, the Earth atmosphere, 
or the Sun light reflected by the Moon. To improve their preci¬ 
sion, those diagnoses can be measured on the averaged line pro¬ 
file which is obtained by computing the weighted cross-correlation 
function (Baranne et al. 1996; Pepe et al. 2002). Note that this tech¬ 
nique conserves the line-profile shape only if the observed spec¬ 
trum is correlated with a binary, box-shaped mask (Eggenberger & 
Udry 2010). The correlation with a non-box-shaped mask such as 
a master-spectrum template (Anglada-Escude & Butler 2012) or a 
synthetic model might smooth out the diagnoses variation, if any, 
and thus decrease the efficiency of the method presented here. 

2.2 Blend simulation 

To test the effect of a contaminating star, we simulated, as described 
in Appendix B, the line profile of a target star with a radial velocity 
of 0 kms^', a FWHM of lOkms^', and a contrast of 30%. Those 
values are typical for a main-sequence star with solar metallicity 
and a slow rotation observed at high-resolution (R > 40000). We 
then blended this main line profile with a contaminating one with 
a FWHM of 70%, 100%, and 130% of the target FWHM, a radial 
velocity shift, 0, from the target star ranging from 0 = Okms^' 
to 0 = 35kms“' and a contrast of 30%. We fixed a flux ratio 
between the target and the blend lines of 1 to 100. Those values 
have been chosen to illustrate the effect of a blend on precise 
spectroscopic diagnosis. Changing their value will not change the 
general behavior presented here, it will just rescale all the values. 
For each pair of radial velocity shift and FWHM of the blended 
star, we measured the radial velocity, contrast, FWHM as well as 
the five line-asymmetry diagnoses presented in Section 2.1 and in 
Appendix A. The results are displayed in Fig. 1. The dash, solid 
and dot-dash lines represent the measured diagnoses as function of 


the radial velocity shift between the two stars for a FWHM of the 
blend of 7 km s“', 10 km s“*, and 13 km s“', respectively. 

We also check the correlation between the measured diag¬ 
noses and the measured radial velocity. The correlation plots are 
displayed in Fig. 1. We discuss the results of these simulations in 
the Sections 2.3 - 2.6. 


2.3 Effect on the radial velocity 


The effect of a contaminating star on the measurement of the radial 
velocity is displayed in Fig. 1. Such perturbation is qualitatively 
similar to the Rossiter-McLaughlin anomaly, where a faint frac¬ 
tion of the line profile is “subtracted” by the eclipsing object. This 
Rossiter-McLaughlin anomaly has been described analytically, e.g. 
in Boue et al. (2013, eq. 15), assuming Gaussian profiles and a 
small eclipsing object (i.e. a planet) compared with the stellar ra¬ 
dius. On the opposite, here, the faint line of a contaminant is added 
to the target one. If we assume that both the target and the contam¬ 
inant have Gaussian profiles and a flux-ratio small enough (i.e. the 
contaminant is much fainter than the target), we can approximate 
the RV-blend anomaly behavior as the opposite of the Rossiter- 
McLaughlin anomaly. Thus, following Boue et al. (2013), we can 
describe the radial velocity-blend anomaly {^rv) by the equation: 


^rv(4>) = 


2crj 


3/2 


0/exp(- 



( 1 ) 


where crj and (T 2 are the Gaussian widths of the two line profiles, 
which depends on the v sin i* and the instrumental resolution, 0 is 
the radial velocity shift between them and / is the flux ratio. 


Using equation 1, we can determine the radial velocity shift 
that maximize the perturbation (4>max) by computing the zero of its 
first derivative: 

d^RV (0) ^ ^ , / 2 , 2 

- — - = 0 => = ± yo"! + 0-2- (2) 

Therefore, a contaminating star affects the measurement 
of the radial velocity in two different ways, depending if |0| is 
lower or greater than \(pmax\- If 101 < l0ma.vL the observed radial 
velocity will be correlated with the radial velocity variation of the 
blending star. On the opposite, if |0| > \^,„ax\, the observed radial 
velocity will be anti-correlated with the radial velocity variation 
of the contaminating star. Note that by definition of the Rayleigh 
criterion, if |0| > |0 „m;c|, the two line profiles are spectroscopically 
resolved and if |0| < \4imax\, they are spectroscopically unresolved. 


Now, if we compute the maximum amplitude of ^rv when 0 = 
lipmaxl in equation 1, we find that: 

(0mflx) = 2 V2 /expl-H. (3) 

<T\ + crl i 2j 

This equation shows that the larger Ci is, due to stellar rotation 
or to lack of instrumental resolution, the larger is the RV-blend 
anomaly signal. 


The RV-blend anomaly displayed in Fig. 1 is similar to the one 
found by Cunha et al. (2013), except that we do not found the sec¬ 
ondary anomaly reported by these authors at large RV-separation, 
which is caused when the contaminant leaves the cross-correlation 
domain. For larger separation, the lines are resolved and this tech¬ 
nique does not provide constraints to planet validation. 
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Figure 1. Top and middle plots: Measured RV, FWHM, contrast, BIS, Vspan, Wspan, BiGauss, Vasy (from left to right and top to bottom) as function of the 
RV offset between the target and the contaminating stars (<j), see text), for a flux ratio in the continuum between the two stars of 1%. Bottom plots: Measured 
con'elation between the radial velocity and FWHM, contrast, BIS, Vspan, Wspan, BiGauss, Vasy (from left to right). The dash, solid and dot-dash lines 
represent the measured diagnoses as function of the radial velocity shift between the two stars for a FWHM of the blend of 7 km s“*, 10 km s“*, and 13 km s“* 
(respectively). 


2.4 EfTect on the FWHM 

The FWHM is affected by the presence of a contaminant. In the 
cases that we simulated here, the variation in FWHM is about two 
times larger than the one in radial velocity (see Fig 1). As for the 
radial velocity, the FWHM varies in case of blended star in two 
different regimes: the FWHM is correlated with ^ for small RV 
separation and anti-correlated for large RV separation. The maxi¬ 
mum of the FWHM-blend anomaly occurs at a different value of 
If the two stars have similar radial velocities {(j) x 0), a small ra¬ 
dial velocity variation of the contaminant will not produce a large 
variation of the FWHM value. The variation is even fainter if the 
contaminant star is rotating fast. 


2.5 Effect on the line contrast 

The contrast of the line profile is also affected by the contaminat¬ 
ing stars. Compared with the other diagnoses, the contrast presents 
only small variations. Moreover, the contrast can be perturbed by 
several instrumental systematics such as residuals of the CCD dark 
current correction, non-linearity of the CCD, sky luminosity varia¬ 
tion caused by the Moon, twilight, or city lights or even sky trans¬ 
parency related to the humidity and the presence of clouds. There¬ 
fore, the contrast of the line profile is not a good diagnosis to con¬ 
strain blended components. The contrast might still be used to high¬ 
light unseen stars in a spectrum (e.g. O-type stars) which would not 
present many stellar lines apart from the Balmer lines, but would 
dilute substantially the spectrum of the contaminant. In this case, 
a narrow line profile with a low contrast could be observed, and 
the low-contrast anomaly might be used to reveal the presence of a 
massive or very fast-rotating star. However, in extremis, this could 


also be mimicked by a metal-poor star. For these reasons, we will 
not mentioned it anymore in the rest of the paper. 

2.6 Effect on the line asymmetry 

It is well known that a contaminating star affects the asymme¬ 
try of the line of the target star if their radial velocity separation 
is small enough. All the asymmetry diagnoses present a similar 
anomaly curve, with different amplitude, except for Vasy which 
present an anomaly curve with an opposite sign (note the reversed 
Y axis in Fig. 1 for the Vasy diagnosis). This opposite shape of 
the Vasy anomaly curve is due to its different definition compared 
with other diagnoses (see appendix A for more details about their 
definition). Compared to the other diagnoses, we observe three dif¬ 
ferent regimes of the asymmetry-anomaly curves presented in Fig. 
1, depending on the FWHM of the contaminating star (FWHM 2 ), 
compared with the target star (FWHM;): 

(i) FWHMi < FWHM 2 : if the contaminating star is rotating 
faster than the target star (see the dot-dash line in Fig. 1), the 
anomaly curves of the BIS, Vspan Wspan, BiGauss and Vasy present 
a shape that is similar with the radial velocity-anomaly curve. The 
main difference is that the maximum anomaly is reached for a 
radial velocity separation of about twice (pmax (see Fig. 1). Such 
separation difference in the maximum of the anomalies is the cause 
of the “figure eight loop” seen in the correlation between the line 
asymmetry and the radial velocity (see Fig. 1). This shape of the 
BIS - RV correlation plot has already been observed by Boisse et 
al. (2011a) for stellar activity, but with an opposite sign (see for 
example their figure 6). 

(ii) FWHM] « FWHM 2 : If the two stars have similar FWHMs 
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(see the solid line in Fig. 1), the asymmetry diagnoses are not very 
sensitives to a contaminant with a small radial velocity separation. 
This configuration could he quite common when the FWHM of 
both stars are dominated by the instrumental resolution and not 
by stellar rotation. A contaminating star with the same FWHM 
than the target star and with a radial velocity variation up to a few 
kms^’ around the target systemic radial velocity will only slightly 
affect the BIS, Vspan, Wspan, BiGauss, and Vasy. In this blind zone, 
the observed radial velocity and asymmetry diagnoses might not 
produce a significant correlation (see also Section 2.9). 

(iii) FWHMi > FWHM 2 : If the target star is rotating faster 
than the contaminating star (see the dash line in Fig. 1), all asym¬ 
metry diagnoses present three regimes: for small values of 0, the 
measured asymmetry is anti-correlated with (j) and thus, with the 
observed radial velocity. Then, for intermediate values of 0, the 
asymmetry diagnoses are correlated with the observed radial ve¬ 
locity. Finally, as in case (i), when the two stars have a large RV 
separation, the asymmetry diagnosis is anti-correlated with (j) but 
correlated with the observed radial velocity. 

2.7 Photon-noise uncertainty 

The various asymmetry diagnoses presented in the previous sec¬ 
tion are not of equal sensitivity to blend scenarios. The diagnosis 
that presents a maximum of sensitivity to a contaminant is Wjp^^, 
then Bigauss, BIS, and finally Vjpa^. This has to be compared with 
the sensitivity to noise in the line profile. For that, we estimated the 
photon noise uncertainty in the various diagnoses. We took a solar 
spectrum observed on the blue sky with a signal to noise (SNR) 
greater than 150 by the HARPS (Mayor et al. 2003) and SOPHIE 
(Perruchot et al. 2008; Bouchy et al. 2009b) spectrographs. For SO¬ 
PHIE, we used two different spectra, one for each instrumental con¬ 
figuration: the high-resolution (HR, R ~ 75000) and high-efiiciency 
(HE, R ~ 39000) modes. For comparison, HARPS has a resolution 
of R ~ 110000. We add white noise to these spectra at the level 
of 10%, 5%, 2%, and 1% of the flux in each pixel, that we gen¬ 
erated 100 times. We ended with 400 spectra for each instrument. 
We cross correlated all these spectra using a G2V numerical mask 
and measured all the diagnoses. We then computed the standard 
deviation of all the diagnoses for each noise level. They all exhibit 
a clear correlation with the radial velocity photon noise. We fitted 
this correlation with a 1-D polynomial and report in the Table 1 the 
value of the slope, for each diagnosis and each instruments. There¬ 
fore, these values should be used to scale the diagnoses uncertainty, 
given the radial velocity photon noise such as: 

CTk = SkO-Rv, (4) 

with CTk the photon noise uncertainty of the diagnosis k, Sk the 
scaling coefficient listed in Table 1, and the radial velocity 
uncertainty. Note however that those values represent only the 
photon noise uncertainty. Those diagnoses are also sensitives to 
instrumental effects (e.g. Boisse et al. 2010; Boisse et al. 2011b; 
Dfaz et al. 2012; Santerne et al. 2014) which might produce 
important additional red noise. Note that the typical photon noise 
reached by SOPHIE and HARPS on stars in the magnitude range 
12 - 16 is presented in Santerne et al. (2011) and compared with 
the expected radial velocity amplitude of the Kepler candidates in 
Santerne et al. (2013b). 

We find that the photon noise uncertainty of the FWHM is 
about a factor of two bigger than the RV, as for BIS and BiGauss. 


Table 1. Photon noise uncertainty on the spectroscopic diagnoses, relative 
to the one of the radial velocity, expressed in kms“*. The contrast is ex¬ 
pressed here in %. 


Spectrograph 

FWHM 

BiGauss 

BIS 

^span 

Wspan 

Vasy 

Contrast 

HARPS 

2.0 

2.1 

2.0 

1.5 

5.8 

5.6 10“* 

11.3 

SOPHIE HR 

2.5 

2.3 

2.3 

1.6 

6.8 

3.6 I0‘“ 

10.2 

SOPHIE HE 

2.5 

2.7 

2.6 

1.8 

8.4 

2.5 10*“ 

7.0 


This agrees with the values usually assumed (Queloz et al. 2009, 
e.g.). As claimed by Boisse et al. (2011a), Vjpan is less noisy than 
BIS by a factor of about 25%. However, the new line-asymmetry 
diagnosis we present in this paper is less robust to noise than BIS, 
its photon noise uncertainty is larger by a factor of about 3. The 
dimension-less diagnosis Vasy present a high-value uncertainty, 
which is difficult to interpret. By comparing the amplitude of the 
Vasy-blend effect in figure 1 with the one of the radial velocity, 
we realised it exhibits a variation AVasy / ARV 3.10^ km^'.s 
larger than the one in radial velocity. However, its photon noise 
uncertainty scaling coefficient, Sv^sy, is at the level of a few 10'° 
km^'.s. Therefore, the diagnosis Vasy seems to be more sensitive 
to the noise than to line-profile asymmetry. Based on that, we do 
not recommend to use Vasy to monitor line-profile asymmetries. 
This is apparently not similar for the contrast, which exhibits in the 
simulation of Section 2.2 a RV-normalised variation at the level of 
Acontrast / ARV = 34%.km“*.s while its photon noise scaling coef¬ 
ficient is at the level of Scontrast = 10%.km“'.s. However, the contrast 
is affected by numerous systematic effects discussed in Section 2.5, 
which make this diagnosis not reliable for the purpose of this paper. 


2.8 Comparison of the diagnoses sensitivity 

We report five different line-asymmetry diagnoses that present 
different sensitivity to blend and to the photon noise. To compare 
them, we normalised the anomaly curves produced in Section 2.2 
with the photon-noise coefficients of the Table 1 . We find that the 
most sensitive diagnoses to line-profile asymmetry are BiGauss 
and Wspan (see Fig. 2). They are followed by BIS and V^pa,,. Being 
a dimension-less diagnosis, Vasy is difficult to compare with the 
other ones. However, as discussed in the previous section, Vasy is 
more sensitive to noise than line-profile asymmetry. Thus, we do 
not consider it further in the end of the paper. Compared to Wspan 
and Vspan, BiGauss does not need a high oversampling of the line 
profile for their fitting procedure to be stable (see Appendix C). 
BiGauss is also the fastest one to compute at a given fit precision. 
Even if we did not compared those diagnoses in the case of stellar 
activity, we expect the same results concerning their sensitivity to 
line-profile asymmetry. For all these reasons, we recommend to 
use BiGauss to monitor line-profile asymmefry. 

We also see in Fig. 2 that the RV-, FWHM-, and asymmetry 
diagnoses exhibit different blend anomaly curves, except when the 
two stars are completely separated spectroscopically. The FWHM 
and the line-profile asymmetry diagnoses do not have a zero slope 
for the same value of (j). Therefore, by combining both the FWHM 
and the line-profile asymmetry diagnoses, one should be able to 
identify unambiguously a blended contaminant, if enough precision 
on those diagnoses is reached. 
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Figure 2. Comparison of the blend curves for all diagnoses but Vasy pre¬ 
sented in this paper. The simulation is the one described in Section 2.2, 
for a contaminating star with a FWHM of 7kms“*. The amplitude of 
the anomaly curves has been divided by the photon-noise coefficients for 
HARPS from Table 1. 


2.9 Illustration with a planet orbiting a companion star 

To illustrate the effect of a contaminating star, we arbitrary 
simulated a binary system composed by a primary star of IMq 
and a secondary star of 0.6Mo in a long-period orbit. In the V 
band, the flux ratio between the two stars is of about 1:30 for 
main-sequence stars with solar metallicity. We assumed a systemic 
RV shift between the two stars of 0o = 1 km s“'. We then assumed 
that the secondary star has a dark companion that produces a RV 
variation with a semi-amplitude K 2 = 50ms“* (this corresponds 
to a Jupiter-mass planet at 200 days or to a Saturn-mass planet at 
5 days). As in the previous section, we simulated the observation 
of such system assuming that the primary star has a FWHMj of 10 
kms“* and the secondary has a FWHM 2 of 7 kms^' (dash line), 
lOkms^' (solid line), and 13kms“* (dot-dash line). For simplicity 
in this example, we considered a circular orbit. The blended radial 
velocities of this system present a nearly circular apparent signal 
with a semi-amplitude of about 1.5 m s“*. The dilution of the signal 
is discussed in Section 3.1. 

The blended signal is nearly but not strictly circular, as 
already pointed out by Wright et al. (2013). As shown in the 
middle panel of Fig. 3, the residuals from a circular fit exhibit 
another signal exactly at half of the initial period. In the present 
simulation, this residual signal has an amplitude of about 0.1% of 
the one of the blended signal (hence 0.15 cm.s“‘). This P/2 residual 
signal should be maximum for 0o = 0 and K 2 = (f/max (eq- 2 - here 
K 2 ~ 6 kms^' which would correspond to a 5Mjupbrown dwarf at 
5 days). In this particular case, the blended radial velocity signal 
exhibits an amplitude of about 100 ms“* and the residual signal 
has a relative amplitude of 6% (hence 6ms“', for FWHMi 2 = 10 
kms“‘). If K 2 > (pmax or if the reflex motion of the secondary star 
does not always satisfy \^\ < along the orbit, then the blended 
signal is clearly not a Keplerian, as already illustrated in Torres et 
al. (2005). Note that, as found by Wright et al. (2013), this residual 
signal could also be at the second harmonic, hence at P/3. 

We also display in Fig. 3 (bottom panels) the correlations be¬ 
tween the relative RV and the relative FWHM, BIS, V^p^^, Wjpa^, 
BiGauss and Vasy. We do not display here the contrast since it 
presents only a variation at the level of ~ 0.3%, which is unde¬ 
tectable with current instrumentation. The blended signal revealed 
a positive-slope correlation between RV and BIS, Vjpan, Wjpan and 





Figure 3. Simulation of a planet orbiting a companion star (see text). Upper 
pane!', blended radial velocities as function of the orbital phase and relative 
to their systemic velocity. Middle panel: relative residuals of the blended 
radial velocities from a circular fit as function of the orbital phase of the 
companion transiting planet. The residuals exhibit a signal at half the ini¬ 
tial period. Lower panels: correlation between the measured RV and the 
FWHM, BIS, Vspan, Wspan, BiGauss and Vasy (from middle to bottom and 
left to right). The different lines represent different FWHM2: 7 km s“* (dash 
line) lOkms^* (solid line), and 13kms“* (dot-dash line). 


BiGauss for FWHM; < FWHM 2 and a negative-slope correlation 
otherwise. If the two stars present a similar FWHM, then the corre¬ 
lation between the aforementioned line-asymmetry diagnoses and 
the RV is very weak, which might be undistinguishable from noise. 
Only the FWHM might present a clear correlation with the RV in 
this case. This shows that it is possible to have a false positive sce¬ 
nario which produce the radial velocity variation compatible with a 
planet but no correlation between the RV and the line-profile asym¬ 
metry diagnoses. It is also possible that a false-positive scenario 
produce a negative-slope correlation between these two quantities, 
which might be misinterpreted as stellar activity. In both cases, the 
FWHM correlates with the RV. Thus, the FWHM is a good diag¬ 
nosis that should be used in complement to other line-asymmetry 
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diagnoses to constrain blend scenarios and infer the planetary na¬ 
ture of a signal. 

3 SENSITIVITY OF RADIAL VELOCITIES TO BLENDS 
3.1 Dilution of the radial velocity signal 

The presence of a blended contaminant dilutes any radial velocity 
variation (Buchhave et al. 2011; Cunha et al. 2013). As exemple, we 
simulated in the previous section a 50 m s^'-signal which was seen, 
after dilution hy another star, as being less than 2 m s“*. To quantify 
this dilution in radial velocity, we evaluated the ratio between the 
amplitude of the blended RV variation (i.e. the one that would he 
observed) and the real one. In the approximation of low flux ratio 
and Gaussian profiles, we can use eq. 1 to estimate the blended 
RV variation, given the RV variation of the contaminating 

star, A0 (i.e. a variation of RV offset). In other words, A^rv is the 
radial velocity variation that would he observed while A0 is the true 
one. We can therefore describe the radial velocity dilution factor 
A^rv'/A<(i as: 



Figure 4. Iso-contours of the RV dilution ld<p as function of the nor¬ 
malised RV shift between the main and secondary stars, and the 

flux ratio between their line profile, /. Solid lines represent positives values 
of d^Rvld<f> while dash lines stand for negative values. Flux ratios greater 
than 0.1 are not displayed here because they do not respect the approx¬ 
imation used in this computation (see text). For this figure, we assumed 
FWHM = (Ti X 2 V2ln(2) = lOkms-'. 
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The figure 4 displays the iso-contours of d^RvId^ as function of 
4’14’max and /. for FWHM = tri x 2 V2 ln(2) = lOkms^'. Note that 
these analytical results are perfectly reproduced by our numerical 
simulation for all / « 1. Three regimes are clearly visible: 


• (j) (titmx- if the two line profiles are unresolved, the RV 
dilution scales as A^RvlAcp oc /. 


• <p X if the two line profiles are shifted exactly by 4'max^ a 
small variation of the contaminating star would not result in a RV 
variation of the blended line profile. In this case, the radial velocity 
is insensitive to faint blend variations. 


• (p > (pmax'- if the two line profiles are resolved, the RV dilution 
scales as A^ry IAtp oc This means that the blended RV variation 

will be anti-correlated with the contaminant variation. If a particu¬ 
lar epoch of the orbit is constrained by other means (e.g. the transit 
epoch is known thanks to photometric data), it would result in a 
blended RV variation in anti-phase (as reported among the CoRoT 
detections, e.g. Cabrera et al. 2009; Erikson et al. 2012; Carone et 
al. 2012; Cavarroc et al. 2012), or with a negative amplitude (e.g. 
Marcy et al. 2014). 

We assume now the opposite case, where the brightest com¬ 
ponent has a RV variation A0' and the faintest component has a 
constant RV in time. In the referential of the brightest star, this RV 
variation will be seen as a RV variation of the secondary star Atp 
such as Atp' = -A0, hence the RV dilution will be: 


A(p' 


A(p 



A0 


0 , 


(V) 


with A^'j^y the blended radial velocity variation given a real varia¬ 
tion of the main star of Atp'. The same approximations are used in 
this equation than for eq. 5. We validated this analytical expression 
with a numerical simulation. Therefore, following the three regimes 


displayed in figure 4, we find that if a planet host is contaminated 
by a secondary star, hence the blended RV variation will be: 

• smaller than the real variation if the two stars are unresolved 

(0 < Ipmax) ; 

• unaffected if the two stars are shifted exactly by (p^ax ; 

• greater than the real variation if the two stars are resolved 

(^ > 0„,a.v). 


An illustration of the first case is the WASP-85 A & B bi¬ 
nary system which hosts a transiting planet (Brown et al. 2014). 
In this system, the unblended HARPS radial velocities present a 
larger amplitude than the blended SOPHIE and CORALIE data. 
Note however that the formalism developed here can not be ap¬ 
plied to this case since the flux ratio between the two stars is not 
low enough. This system could have been resolved though numer¬ 
ical simulation, as done in the FASTIS software, but the expected 
variation of BiGauss and FWHM is at the level of 1.5 ms“* and 
10 m s“' for SOPHIE (respectively), which is much below the pre¬ 
cision on this star. Other illustrations of radial velocity dilution by 
an unresolved contaminant is Kepler-14 (Buchhave et al. 2011) and 
KOI-1257 (Santerne et al. 2014). As discussed in the later paper, 
since nearly half the stars are in binary system, there might be much 
more diluted radial velocity cases than we have found so far. 


3.2 Radial velocity sensitivity to physical companions 

In the previous section, we presented the dilution of a given 
radial velocity signal as function of the flux ratio and the RV shift 
between the blended stars. To interpret this radial velocity dilution 
as a sensitivity to blends, we related the flux ratio to the physical 
properties of the blended stars. In that case, we tested the sensitiv¬ 
ity of radial velocities to a companion star, i.e. bounded with the 
primary star. For this scenario, we considered a primary star with a 
mass of Ml = IMq and a secondary one with M 2 6 [0.1; 1.0] Mq. 
We set that both stars have solar properties ([Fe/H] = 0 dex, age 
= 5Gyr), the same distance, a r/sinf* of 3kms“*, and a systemic 
radial velocity shift of 00 = 4 km s“*. We simulated a faint radial 
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velocity variation of the secondary star, caused by the presence of 
a sub-stellar, non-emitting, companion. The RV variation (Ai^) of 
this secondary star is kept relatively small (we assumed a circular 
1-Mjup companion at 365 days of period, hence K 2 < 130 ms“*). 

We tested this scenario for the SOPHIE spectrograph (Perru- 
chot et al. 2008; Bouchy et al. 2009b), for which the line profile has 
been calibrated by Boisse et al. (2010). We used these calibrations 
(also reported in Appendix B) to simulate this system for the two 
instrumental configurations: the High-Efficiency mode (HE) with a 
spectral resolution of ~ 39000 and the High-Resolution mode (HR) 
with a spectral resolution of ~ 75000. We also considered the two 
masks G2 and K5 that are used to compute the cross-correlation 
function in the online pipeline. For each M 2 , instrumental mode 
and mask, we simulated the CCFs of both stars, add them, and 
measure the blended RV variation that we divided by A0. 
For that, we used the Dartmouth evolutionary tracks (Hotter et al. 
2008) to estimate the stellar parameters (radius, Tefr, luminosity) 
and the BT-SETTL (Allard et al. 2012) spectral atmosphere models 
to estimate the flux of both stars within the spectrograph bandpass 
(assumed to be similar to Johnson-V) and their (B-V) color index 
needed to derive their CCF parameters (see Appendix B). 

We display in figure 5 (upper panel) the ratio A^rv/A 0 as 
function of the mass of the bounded secondary star, host of the 
sub-stellar companion, for both instrumental modes and correla¬ 
tion masks. In this configuration and with a 1 ms“* precision in 
the radial velocity we can constrain the signature of a blended 
K 2 = 100 ms“' RV signal on a companion with a mass down 
to ~ 0.5Mo. To detect the same signal on a O. 3 M 0 , one needs a 
precision at the level of a few lOcms^'. As expected, the blended 
RV variation decreases as the secondary star become smaller and 
thus fainter in the V-band (the flux ratio decreases). We observe 
that for M 2 « Ml, the radial velocity signal is larger in the mask 
K5 than in the mask G2. This is the so-called mask effect already 
reported in Bouchy et al. (2009a). The K5 mask has more lines 
than the G2 one. Thus, the CCF of a blended G2-I-K5 spectrum will 
present a slightly larger flux ratio if correlated using all the lines of 
a K5 star compared with the lines of a G2 star. Since a larger flux 
ratio implies a larger blended signal, the K5 mask is more sensitive 
to K and M dwarfs than the G2 mask. The difference of radial 
velocity as derived by the two masks therefore provide constraints 
on the mass of the contaminant star compared with the one of the 
primary star. 

In the figure 5, we also observe that the HE mode of SOPHIE 
is more sensitive to blends than the HR mode. In eq. 6 , the factor 
d^Rvld(/> scales as 1 - {4>l4>maxf for (j) < cpmax (unresolved compo¬ 
nents). Since the instrumental resolution of the HE mode is lower 
than the HR one, the Gaussian widths of the stars are larger in HE 
than in HR, which makes (pmax larger, and thus d^RvIdij) is also 
larger. Therefore, decreasing the instrumental resolution increases 
asymptotically the sensitivity of radial velocities to unresolved 
stellar components. On the other hand, for (j) » (ji^ax (i-e- resolved 
stellar components), the factor d^RvId^ is dominated by the term 
exp(-l/2(0/(^,„„.t)“). In that case, increasing the instrumental 
resolution also increases the sensitivity of radial velocities to 
resolved stellar components. 

To constrain the presence of a stellar contaminant in radial 
velocity data, a good test is either to use spectrographs with 
different instrumental resolution or to artificially decrease the 




Figure 5. Relative amplitude of a blended radial velocity variation, 
A^rv / A^, as function of the mass of the blend host M2 for the considered 
scenario (see text). Upper panel: The four lines display this relative ampli¬ 
tude for both the SOPHIE HE (black lines) and HR modes (grey lines) using 
cross-correlation mask of a G2 (solid lines) and K5 dwarf (dash lines). This 
assumes cpa = 4kms“'. Lower panel: Iso-contours of this relative ampli¬ 
tude as function of the blend mass M2 and the distance modulus difference 
ixd (see text), for the SOPHIE HE mode and the G2 mask. The hatched 
region represents contaminant stars brighter than than the target in the V- 
band. 


instrumental resolution by convolving the spectra with a Gaussian 
profile of increasing width. Then, if the amplitude of the radial 
velocity signal varies significantly with the instrumental resolution, 
the line-profile is thus blended. Note that the opposite is not always 
true: since this is an asymptotic effect, no amplitude variation as 
a function of the instrumental resolution does not imply that the 
line profile is unblended. At the most extreme case, when the two 
stellar components have perfectly the same radial velocity (0 = 0 ), 
the factor d^Rvldc/) scales as So, in that case, if 

both Cl and cr 2 are dominated by the instmmental resolution and 
not by the stellar rotation or if both stars have the same u sin i* , 
the factor d^Rvldc/) is independent of the instmmental resolution. 
However, if only one star is rotating fast compared with the 
instmmental resolution, this test should be conclusive, even if the 
stars have a similar systemic radial velocity. 

Note that high-precision spectrographs in the near infrared, 
like SPIRou (Delfosse et al. 2013) or CARMENES (Quirrenbach 
et al. 2014), should also provide further constrain on the presence 
of low-mass stars, by improving the flux ratio between the primary 
and the lower-mass companion stars. 
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3.3 Radial velocity sensitivity to background companions 

A similar test can be done by relaxing the assumption that the pri¬ 
mary and contaminating stars are bounded. This scenario has to as¬ 
sume that both stars are chance-aligned in both the plane of the sky 
and in the RV space, which is quite unlikely (Cunha et al. 2013). 
Even though, we can quantify the constraints that the radial ve¬ 
locities could provide. For that we repeat the previous simulation 
by allowing the distance ratio between the two stars to vary in the 
range Ad e [-6; -l-6], with Ad the distance modulus difference as 
defined in Fressin et al. (2011): 

Ad = 51og,o(^j, (8) 

where Di and D 2 are the distance of the primary and contaminating 
stars, respectively. For this test, we set 0o to zero for simplicity, 
which gives optimistic values of the radial velocity constraints. 
The lower panel of figure 5 shows fhe iso-contours of A^«v/A0 as 
function of the mass of the host star (i.e. the contaminating star) 
and the distance modulus difference Ad for fhe G2 mask. Since 
00 = 0, bofh insfrumenfal modes give the same constraints here. In 
this configuration, with a 1 ms“' precision on SOPHIE, one could 
constrain a perfectly chance-aligned system with K 2 = 100 ms“’ 
on a IMq host with Ad < 5 or a O.VMq host with Ad < 2. This 
map of the radial velocity sensitivity to blends could be compared 
with the constraints on background systems provided, e.g. by 
transit photometry, if any. From the photometric constraints on the 
possible false positive parameters, one could estimate the expected 
amplitude of the radial velocity signal of such false positive using 
this technique. Then, it would be possible to optimise the radial 
velocity follow-up of transit candidates as function of the precision 
of the available spectrographs and their expected efficiency to 
further constrain, or reject, a given false-positive scenario. This 
could permit to minimize the needs of radial velocity spectrograph 
to screen out false positives of the future transit space missions 
TESS (Ricker et al. 2014) and PLATO (Rauer et al. 2014). 


4 BENCHMARKING THE BLEND MODEL: THE CASE 
OE HD16702, A NEARLY FACE-ON BINARY 

To validate our model of blended spectroscopic diagnosis, we 
tested it on one case which has other, independent constraints. 
This case is the system HD16702 revealed in Dfaz et al. (2012). 
Observed by SOPHIE, this system exhibit a radial velocity signal 
compatible with a brown dwarf companion. A careful analysis of 
the line bisector only, with a preliminary version of this model, 
revealed that the system is a nearly face-on binary. The same 
authors also reported the marginal detection of an astrometric 
signature in the HIPPARCOS data compatible with this scenario. 
Recently, Kolbl et al. (2015) developed a new technique to detect 
faint components in the Keck/HIRES spectra. They applied 
this new technique to HD16702 and confirmed the bisector and 
astrometric results. In this section, we reanalyse the SOPHIE 
data of HD16702 with our blend model, by fitting simultaneously 
the radial velocities, FWHM, and BiGauss derived in the two 
correlation masks G2 and K5. 

We assumed a scenario of a nearly face-on binary. We fitted 
all the data simultaneously using the FASTIS software (see Dtaz et 
al. 2014, and reference therein) which uses a Markov Chain Monte 
Carlo (MCMC) method. We used the Dartmouth evolutionary 


tracks (Dotter et al. 2008) to estimate the stellar parameters and the 
BT-SETTL models (Allard et al. 2012) for the stellar atmospheric 
spectrum. The exhaustive list of parameters used for this analysis 
is reported in Table Dl. We ran 15 chains randomly started from 
the joint prior distribution and all chains converged toward the 
same solution. We thinned and merged them to derive the posterior 
distribution. The derived value for the parameters and their 68.3% 
uncertainty are also reported in Table Dl. 

We display in figure 6 the phase-folded data considered in this 
analysis together with the best-fit model. While the radial veloc¬ 
ity and BiGauss present a significant variation, the FWHM does 
not present a clear variation. This absence of significant variation 
is still used in the fitting procedure to constrain the v sin t* of the 
secondary star. In Dtaz et al. (2012), only the bisector was anal¬ 
ysed and it was not possible to break the degeneracy between the 
FWHM of the secondary and its mass (see the fig. 9 in Dfaz et 
al. 2012). In this analysis, our model was able to constrain the 
V sin I* of the secondary star, as displayed in the figure 7. This was 
possible thanks to the additional constraints from the FWHM. We 
also ran the same analysis but for each correlation mask indepen¬ 
dently and found similar results, with slightly larger uncertainties 
as expected. Since the degeneracy between the mass of the sec¬ 
ondary star and its usint* is broken, we can further constraints 
the mass to a value of M 2 = 0.35 ± 0.03Mo. Note that this un¬ 
certainty is only the statistical one. Uncertainties from the stellar 
models or from the blend model related to the calibration of the 
CCF are not taken into account here. This value is fully compatible 
with the values found by Dtaz et al. (2012) and Kolbl et al. (2015) 
of M 2 = 0.40 ± 0.05Me and M 2 ~ O. 3 M 0 (respectively). This 
mass is also fully compatible with the one derived by the HIPPAR¬ 
COS astrometry of M 2 = 0.55 ± 0.14Mo, which was used as prior 
for this analysis (see Table Dl). Astrometric measurements from 
the GAIA space telescope (Prusti 2012) will permit us to confirm 
this result. This allows us to refine the inclination of the binary to 
iorb = 9.4 ± 0.5°. The RV-blend model we developed in the frame¬ 
work of the FASTIS validation tool is therefore able to reproduce 
perfectly the spectroscopic data of this blended, nearly face-on bi¬ 
nary. 


5 CAVEATS 

We would like to point out a few caveats in our model and results 
presented in this paper: 

• The calibration of the SOPHIE CCF that we are using is 
correct for (B-V) in the range of [0.43; 0.98] while we are using 
it also for stars outside this domain. We assumed that it could 
be extrapolated without a significant impact on our results. We 
think this is not a limitation of our model because the parameters 
of the simulated line profile have only a second-order effect 
compared with the flux ratio. A new calibration of the SOPHIE 
CCFs with a wider range of (B-V) is mandatory to overcome this 
caveat. Note that calibrating the CCFs with masks corresponding 
to F and M dwarfs would improve the mask effect and thus 
the capability of this model to constrain higher or lower mass 
stars, respectively. Calibrating the CCF of other spectrographs 
in the visible but with diiferent spectral resolution, like HARPS, 
HARPS-N, CORALIE, or in the near infrared could also improve 
the constraints. Note also that at a given (B-V) and [Fe/H], 
the surface gravity should also change the parameters of the 
CCF, which is not accounted for in the current calibration. The 
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Figure 6. Results from the analysis of SOPHIE data of the nearly face-on 
system HD16702. From top to bottom: phase-folded radial velocity veloc¬ 
ities, FWHM, and BiGauss of this system, observed by the SOPHIE in the 
HR mode and derived using a G2 (black points) and K5 (open points) cor¬ 
relation mask. The data are plotted together with the best model found in 
the PASTIS analysis and their residuals. For clarity, an arbitrary offset has 
been set to the K5 data and models in the three plots. 


CCFs of evolved stars should also be considered for the calibration. 

• When we analyse the radial velocity, FWHM and line- 
asymmetry diagnosis from different correlation masks simulta¬ 
neously within the PASTIS software, we compute the likelihood 
for each individual dataset and multiply them to get the marginal 
likelihood. This assumes that the dataset are independent from 
each other, which might not be the case. We assume however that 
this independence of the dataset, if any, does not affect significantly 
our results. 

• We model the line profile of the stars with a Gaussian pro¬ 
file. However, as soon as the stars are rotating faster than the in¬ 
strumental resolution, their line profile differs significanfly from a 
Gaussian profile. Our approximation might lead to over- or under¬ 



Secondary mass [M,.] 

Figure 7. Posterior distribution of the v sin i* of the secondary star in the 
HD16702 system as function of its mass. The regions display the 68.3%, 
95.5% and 99.7% confidence intervals (from dark to light greys) on these 
parameters. 


estimation of the spectroscopic diagnoses. Nevertheless, even if the 
values significantly change, we expect that the overall behaviour of 
the spectroscopic diagnoses is not affected by stellar rotation. 


6 CONCLUSION AND DISCUSSION 

In this paper, we present the radial velocity model developed to 
constrain exoplanet blend scenarios with the PASTIS software. 
Using this model, we investigated how a blended stellar contam¬ 
inant can perturb the measurement of the various spectroscopic 
diagnoses of a target star. We estimated the photon noise precision 
of these spectroscopic diagnoses which allowed us to conclude 
that BiGauss, as defined by Nardetto et al. (2006), is the most 
sensitive diagnosis to line-profile asymmetry. By exploring the 
parameters of both the target and contaminating star, we find 
that a contaminant might produce different FWHM and bisector 
variations but similar radial velocity signals. In particular, if both 
stars have the same FWHM, their blended line profile does not 
exhibit strong asymmetry and can not be monitored with the 
bisector-like diagnoses. However, in that case, the FWHM should 
still present large variations that could be detected with stabilised 
high-resolution spectrographs. Therefore, the absence of bisector 
variation alone can not be used to firmly establish the nature of 
a planet detected by radial velocity. Both precise bisectors and 
FWHMs should not exhibit variation to exclude most astrophysical 
false-positive scenarios. 

More interestingly, we find that if the FWHM of the con¬ 
taminating star is smaller than the target one, the blended line 
bisectors will be anti-correlated with the radial velocities. This 
anti-correlation is usually produced by and interpreted as spots on 
the stellar surface (Queloz et al. 2001; Boisse et al. 2011a, 2012; 
Dumusque et al. 2014). Therefore, a bona-fide planet in a blended 
binary could reproduce the signal in the RV, FWHM and bisector 
of an active star. Such system would produce a hue-negative 
scenario, which might still be identified using other spectroscopic 
diagnoses like the variation of the S -index or of the logfR^,^). 

Our simulation also confirms the results of Wright et al. 
(2013), that the blended signal of a circular radial velocity varia¬ 
tion is not strictly circular. The residuals of a blended circular orbit 
exhibit another signal at the first or second harmonic (P/2 or P/3, 
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respectively) with a relative amplitude up to a few percent of the 
initial blended signal. 

In Section 3, we developed a new formalism to estimate 
the amplitude of the radial velocity dilution caused by a blended 
(brighter or fainter) star. It shows that, for example, a contaminat¬ 
ing eclipsing binary spectroscopically unresolved with the target 
star (i.e. their RV shift is smaller than their joint profile widths) 
might mimic a radial velocity signal in phase with the eclipse 
ephemeris. On the other hand, if the two stars are resolved (i.e. 
their RV shift is greater than their joint profile widths), the blended 
radial velocity signal would be in anti-phase with the eclipse 
ephemeris. This formalism also revealed that if a radial velocity 
variation of a target star is blended with another fainter component, 
the later one would actually boost or dilute the amplitude of the 
signal, depending if the two stars are resolved or not (respectively). 

Using this formalism, we also discussed the sensitivity of 
the SOPHIE spectrograph to detect physical and background 
companions to the target star. By comparing the two instrumental 
configurations of the SOPHIE spectrograph, we realised that the 
instrumental configuration with the lowest spectral resolution was 
more sensitive to unresolved contaminants than the configuration 
with the highest spectral resolution. This is however the opposite 
for a resolved contaminant. Therefore, by using different spec¬ 
trographs with different spectral resolutions, one might constrain 
the presence of a contaminating star. Comparing the amplitude 
of signals from spectrographs with different spectral resolution 
provides another test to detect and constrain the presence of a 
stellar contaminant blended in the line profile of the target star. 
This can be done also by convolving the observed high-resolution 
spectra with a Gaussian profile of increasing widths. 

The method described in this paper makes use of the 
spectroscopic diagnoses measured precisely on high-resolution 
spectra. For that, the spectrograph should reproduced as closely 
as possible the line profile. This is the case of stabilised, fiber-fed 
spectrographs like SOPHIE and HARPS. The various diagnoses 
could be measured simultaneously with the radial velocity with no 
extra observational cost. The detection of line-profile variations at 
the early stage of the observation could actually avoid wasting tele¬ 
scope time observing false positives rather than planets. However, 
with slit-fed spectrographs, such as HiReS@Keck, the shape of the 
line-profile is correlated with the slit illumination profile, which 
strongly limits the measurement of precision spectroscopic diag¬ 
noses (Wright et al. 2013). The use of the iodine cell allows one 
to monitor the instrumental profile needed to measure precision 
radial velocities (Butler et al. 1996), but limits the determination of 
the line profile shape. In such case, applying the method described 
here is impossible, and other observations are needed to firmly rule 
out blend scenarios, which has an extra observational cost. For 
these reasons, stabilised fiber-fed spectrographs are more suitable 
for radial velocity surveys and the follow-up of transiting planet 
candidates. 

It was already known that high-resolution spectroscopy 
can constrain the presence of resolved systems down to a few 
magnitude fainter than the target star (e.g. Borucki et al. 2012). 
We show in this paper, that it is also the case for spectroscopically 
unresolved systems. However, in the latter case, one should 
however analyse carefully the line profile variation. In the con¬ 
text of validation of small transiting exoplanets, this provides 


complementary contraints to photometry, which cannot exclude 
bright contaminants. Some example of constraints provided by 
high-precision transit photometry can be found in, e.g. Torres et al. 
(2011, 2015); Fressin et al. (2011, 2012a,b); Borucki et al. (2012, 
2013). In all these cases photometric constraints can not reject the 
presence of a planet transiting a background, chance-aligned star 
brighter than Ad x 6, which is the domain where spectroscopy is 
the most sensitive. In the scenario of a background system, it is 
quite unlikely that such system would be chance-aligned in both 
the plane of the sky and the RV space. 

This is not the case for a contaminant physically bounded 
with the target star, which has a high probability of having a 
small radial velocity shift. Since the physical companion is most 
likely blended with the target star, searching for a second set 
of resolved lines in the spectra, as done in the aforementioned 
papers, is not expected to provide strong constraints to this kind of 
false-positive scenario. Therefore, without a more careful check of 
their line-profile variation, as presented in this paper, some of the 
planets validated in these papers might actually be larger planets 
transiting a star physically-bound with the target. This is still true 
even if a significant radial velocity variation is detected (Santerne 
et al. 2014). 

The spectroscopic diagnoses we presented here and the 
technique developed by Kolbl et al. (2015) can constrain false 
positive scenarios which are blended in both the plane of the sky 
and the radial velocity space. Another way to constrain blend 
contaminant would consist in searching for abnormal stellar lines 
that do not correspond to the spectral type of the target star. Indeed, 
if a blended contaminant has a different effective temperature than 
the target star, there might be some lines (e.g. molecular lines) 
that are only visible in the spectrum of the coolest component. By 
checking for extra lines in the spectrum of a target star, one could 
also constrain the presence of a contaminating star. This can also 
be done by cross-correlating the observed spectrum with the mask 
of a M dwarf for which all the lines that are common in both a 
M-dwarf and G-dwarf masks have been removed. By doing that, 
if the difference of Teff is large enough, the resulting CCF should 
only present the contribution from the contaminating star. 

As a benchmark test for the blend model we used in this paper, 
we analyse with FASTIS the SOPHIE data of the HD16702 system. 
This system is a blended nearly face-on binary which exhibit a large 
radial velocity and bisector variations, but no significant FWHM 
variation. Our model was able to reproduce all the data and to pro¬ 
vide constraints on the secondary star that are compatible with the 
ones derived with the HIPPARCOS astrometry (Diaz et al. 2012) 
and the technique of Kolbl et al. (2015). 
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APPENDIX A: T HE SPECTROSCOPIC DIAGNOSES 


AI Radial velocity, EWHM and contrast 

The line profile ('P) of a dwarf star with a low v sin i* can be de¬ 
scribed by a Gaussian profile (at the spectral resolution of RV- 
dedicated spectrographs, i.e. between 40 000 and 120 000): 


S (x) = Co 


1 - c exp 


1 (x - Xof \ 

2 0-2 j 


(Al) 


where x is the radial velocity, Cq is the flux in the continuum, c 
is the contrast of the line, xq is the radial velocity of the line, and 
cr represents the width of the line. For fast-rotating stars, the line 
profile can be modeled following the procedure described in San- 
teme et al. (2012a) and references therein. For Gaussian profiles, 
the FWHM is related to cr by: 


FWHM = o- X 2 V21n2 ^crx 2.3548 (A2) 


By fitting the line profile by a Gaussian function, one can there¬ 
fore measure the radial velocity, the contrast and the FWHM of the 
stellar line profile. 
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Figure Al. (Upperpanel) Synthetic line profile (thick and thin black lines) 
contaminated by a diluted star (dotted-red line) with a line contrast of 10% 
of the target star, equal FWHM and shifted by lOkms^*. The green line 
displays the Gaussian fit (eq. Al) to the blended line profile, the magenta 
line displays the gaussian fit to the top of the line profile (eq. A4) while the 
blue line displays the Gaussian fit to the bottom of the line profile (eq. A5). 
The solid red line displays the asymmetric Gaussian fit as defined in eq. 
A7. Cyan and yellow lines display the bisector within the limit defined to 
compute the bottom and the top parts (respectively) of the BIS. The dashed 
blue and red lines display the Gaussian fit to the blue (eq. A12) and red 
(eq. AI3) wings of the line profile. The vertical grey dashed and dot-dashed 
lines represent the ±1(T and ±3cr of the line profile, as measured by the 
symmetric Gaussian fit. (Lower panel) Residuals of the fit of the synthetic 
blended line profile. Displayed lines correspond to the residuals of the fits 
defined for the upper panel. 


A2 BIS 

To monitor the asymmetry of the line profile several diagnoses have 
been defined. Toner & Gray (1988) defined the bisector as “the 
midpoints of horizontal line segments bounded by the sides of the 
line profile”. This series of midpoints is not convenient to interpret 
numerically for each observation. Some diagnoses were therefore 
suggested in order to describe the asymmetry of the line profile in a 
single value. This was the case of the bisector curvature, suggested 
by Gray (1997) which consists in computing the “velocity span of 
the top portion minus the velocity span of the bottom portion of the 
bisector” (Gray & Hatzes 1997). Since this method considers only 
one measurement in the top and bottom parts of the bisector, it is 
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very sensitive to the noise present in the line profile. To improve 
this line asymmetry diagnosis, Queloz et al. (2001) used the Bisec¬ 
tor Inverse Slope (BIS) defined as the velocity span between the 
average of the top and bottom parts of the bisector: 

BIS = BIS,„p - BlSbottom, (A3) 

where BIStop is the average of the bisector between 60% and 90% 
of the total contrast of the line profile and BlSbottom is the average 
of the bisector between 10% and 40% of the total contrast of the 
line profile, where 100% corresponds to the continuum and 0% to 
the minimum of the line profile (see Fig. Al). We find that sev¬ 
eral different limits have been used in previous studies for the top 
part of the BIS: 55% - 85% according to Queloz et al. (2001), 55% 
- 90% according to Dali et al. (2006) and 60% - 90% according 
to Figueira et al. (2013). Since different BIS definitions will lead 
to different results, a unique definition of how to compute the BIS 
should be defined to avoid wrong interpretation when comparing 
BIS data. We propose here to keep the values 60% - 90% as defini¬ 
tive values for BIS computation since they are symmetric with the 
bottom values around the half maximum of the line. 


A3 V 

’ span 

To improve the precision of the bisector, Boisse et al. (2011a) sug¬ 
gested the Vspan diagnosis. Instead of computing the average mid¬ 
point of the horizontal line segment which can be noisy for low- 
SNR spectra, the idea is to compute the radial velocity difference 
measured by fitting two Gaussian functions, one to the top and one 
to the bottom parts of the line profile. The Gaussian fit is expected 
here to be less sensitive to noise in the stellar line than the tradi¬ 
tional bisector. The limits between the top and the bottom of the 
line profile are defined as the ±lcr-limit from the measured radial 
velocity (xq). The two sub-line profiles (T'top and T'bottom) are thus 
defined as (see Fig. Al): 

= T([-oo,xo-lo-]U[xo-Ho-,-too]) (A4) 

'^bottom = 'P ([-oo, Xo - 3cr] U [Xo - Icr, Xq + Icr] 

U [xo -I- 3a-, -l-co]) (A5) 

The diagnosis Vspj^, is therefore defined as: 

Vspan - I (A6) 

where xo,„p and are the radial velocities measured to the top 

of the line profile ('Ttop> eq. A4) and to the bottom of the line profile 
('Fbottom, eq. A5), respectively. 


A4 BiGauss 


Nardetto et al. (2006) defined another diagnosis (called BiGauss) 
to monitor line asymmetry of pulsating stars. It consists in fitting 
an asymmetric Gaussian function to the line profile (see Fig. Al): 


^asy (^) “ 
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where x < Xq 
where x > Xp. 


(A7) 


where A is the percentage of asymmetry (relatively to cr). The 
BiGauss diagnosis is defined as the difference between the radial 
velocity measured with a symmetric Gaussian function (eq. Al) 
and the asymmetric BiGaussian function (eq. A7), as a proxy of 
the line asymmetry: 


(A8) 


A3 V 

T asy 


A fourth diagnosis was defined recently by Figueira et al. (2013) to 
monitor asymmetry in line profiles caused by stellar activity. This 
diagnosis, called Vasy, consists in computing the differential spec¬ 
tral information, as defined by Bouchy et al. (2001), between the 
red and the blue wings of the line profile: 


V = 

'' asy 


'Li{w7 


' W‘ 


hlue\ 


X Wi 


2, (V, 


(A9) 


where W, are the weight of the point computed at the flux level i in 
the line profile (eq. 8 of Bouchy et al. 2001): 


W,- = 



% + <tI 


(AlO) 


where T, and T, are the wavelength and flux of the spectra, respec¬ 
tively, and o-p the detector readout noise. and VV'"'* are the 
blue and red line profile wing information (respectively), and IT, 
their average. If computed in the radial velocity space on the CCF, 
the equation AlO needs to be changed to (Figueira et al., in prep.): 




^2( SCCFi f 
\ dRVi ) 

CCFi -I- o-^ ’ 


(All) 


where RVi and CCFj are the radial velocity and the flux in the CCF, 
respectively, and c is the speed of light. 


A6 3V span 

Four methods have been developed so far to estimate the asymme¬ 
try of a spectral line, because it’s a key diagnosis to correct radial 
velocities from stellar activity (Boisse et al. 2011a), from blended 
contaminants (Santos et al. 2002; Diaz et al. 2012), and from instru¬ 
mental effects (Boisse et al. 2010). Therefore, more sensitive is the 
diagnosis to asymmetry, better the correction or the constraints are 
expected. For that propose, we define here a fifth diagnosis to mon¬ 
itor line asymmetry. It consists on fitting two Gaussian functions to 
the blue and red wings of the line profile (see Fig. Al): 

‘Pbiue = 'I'([-cx.,Xo]U[xo-b3o-,-boo]) (A12) 

= 'I'([-cx.,Xo-3rr]U[xo,-boo]) (A13) 

We can thus defined the asymmetry measurement of this fifth diag¬ 
nosis, that we call Wspa,,: 

Wspan = (A14) 

where xq^^^^ and xq,.^^ are the radial velocity derived from the fit of 
the blue (T'biuei eq. A12) and red (T'redi eq. A13) wings of the line 
profile, respectively. 

In Fig. Al, we illustrate these five diagnoses by modeling 
the line profile of a star, with a FWHM of 10 km s“' and a con¬ 
trast of 30%. We contaminated it with a star with a line FWHM 
of lOkms^', a contrast of 3%, and shifted in radial velocity by 
-blOkms^' from the main star. This example was chosen to present 
a huge line-asymmetry, that can be seen by eye. 


APPENDIX B: CCF MODELING 

To simulate CCFs, we use the equations (B.l), (B.2), (B.3), (B.4) 
and (B.5) from Boisse et al. (2010) for SOPHIE. We first compute 


BiGauss — XqU^q “ -“^oU^o 
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Table Bl. Coefficients of equations B1 & B2 

instru. mode a fS y 6 A 

SOPHIE-HE 10.52 22.56 22.37 6.95 1.64 

SOPHIE-HR 9.90 22.56 22.37 6.95 1.73 

Table B2. Coefficients of equations B4. 

SOPHIE 

as bs cs 

G2 -0.9440 1.4992 3.8807 

K5 0.2615 2.2606 3.9553 


the FWHM of the Gaussian function from the v sin i* of the star 
following the equation: 


/ V sin i* ^ + AVn _ r—— 

FWHM=d - ^5 -- •2V21n(2) (Bl) 

where A is a scaling factor listed in Table Bl. cto is related 
to the color index (B - V) in the range 0.4 - 1.3, similar to FGK 
dwarfs: 

a-o = a-p(B-V) + y{B - Vf - 6{B - V)^ (B2) 

where the coefficient a, p, y, and S are given in Table B1 . We then 
compute the expected contrast C of the CCF : 


APPENDIX C: COMPUTATION TIME AND NUMERICAL 
PRECISION AND ACCURACY OF THE MODELS 

Measuring the diagnoses to a large sample of observations, or 
through a large MCMC analysis takes time. The computation 
speed as well as the precision and accuracy of the fitting procedure 
should be part of the equation when defining the best diagnosis 
to use. These elements should also depend on the number of data 
points sampling the line profile. To quantify them, we simulated a 
blended Gaussian line profile. We defined a target line profile with 
a FWHM of 10 kms^*, a contrast of 30% and a contaminating 
line profile with a FWHM of 12 kms^*, a contrast of 30%, a 
flux ratio of 1/100 compared with the target star, and a constant 
radial velocity shift of lOkms^'. We used different sampling 
rate for the simulation of this blended line profile, ranging from 
lOcms^’ to 1 kms“*. Then, for each sampling size, we simulated 
100 times the same system, by shifting the blended line profile 
randomly within a uniform distribution with a width of 1 kms“*. 
By doing that, the blended line profile was sampled differently in 
all the simulations. On these synthetic line profile we measured 
the various spectroscopic diagnoses aforementioned as well as the 
time taken to compute them using the optimize.leastsq module 
embedded in the scipy package in python2.7. We used the default 
values for this function, except for the xtol parameter that we 
sed to 10“^. This experiment was done with an Intel Xeon CPU 
cadenced at 2.27GHz. Note that the evaluation of the symmetric 
and asymmetric Gaussian functions (used to derive the radial 
velocity, the FWHM, BiGauss, Vspam and W^pan) is implemented 
within FASTIS in a C-l-l- module, together with their Jacobian 
functions, in order to speed up the fitting procedure. 


'W idlWT) 

' FWHM ’ 


(B3) 


where "W is the area of the CCF, expressed as function of the metal- 
licity [Fe/H] and the color index (B —V) \ 


log HU 


[FelH]+as-bs{B-V) 

Cs 


(B4) 


where the coefficient as, bs, Cs are listed in Table B2. These 
relations were calibrated for 0.43 < {B - V) < 0.98 and for 
-0.47 < [Fe/H] < 0.44 for SOPHIE. 


The line profile are assumed to have a Gaussian profile. We 
therefore compute the average line profile for each star considered 
(CCF,) as follow. 


CCF,- (x) = Ti 


1 - C, • exp 

C 

1 

1 

/ FWHMi U 


'v V2V21n(2)J li 


(B5) 


where F) is the flux of star in the continuum integrated within 
the spectrograph bandpass that we assumed to be equivalent to the 
Johnson-V band for SOPHIE. The blended synthetic CCF is there¬ 
fore: 


CCF(x) = 


I,iCCF.(x) 

Z.T 


(B6) 


We then use this synthetic CCF to measure the various spec¬ 
troscopic diagnoses by fitting it with a Levenberg - Marquardt al¬ 
gorithm as implemented into the Scipy package (optimize.leastsq). 


Figure Cl displays the mean time to compute the various di¬ 
agnoses as function of the sampling size. Note that this simulation 
has been performed with a line profile sampled from -50kms“' 
to -f50kms“' around the target star. Decreasing or increasing the 
wavelength or radial velocity range is expected to scale down or 
up (respectively) these results. As expected, the computation time 
increase exponentially when decreasing the sampling size. The 
diagnoses BIS and Vasy which does not require any minimisation 
algorithms are the fastest ones. The diagnoses Vspan and Wspan 
which require two fitting procedure, take twice more time than 
the measurement of the radial velocity, FWHM, and BiGauss, as 
expected. 

The precision and accuracy of the fitting model should also 
depend on the number of points in the line profile. Thus, we also 
investigated the relative precision and accuracy of our models as 
function of the sampling size. We estimated the precision of the 
diagnoses by computing the RMS of the measurements for the 100 
simulations aforementioned. The accuracy has been estimated by 
comparing the median value of the 100 simulations at different 
sampling size with the one with a sampling size of lOcms^*. 
We therefore assumed here that the lOcms^* sampling has a 
negligible noise, compared with the other simulations. We then 
normalised the precision and accuracy with the median value 
of the asymmetry found for a sampling size of lOcms^*, to 
derive the relative precision and accuracy. Figure Cl displays the 
relative precision and accuracy to measure the various asymmetry 
diagnoses as function of the sampling of the line profile. We find 
that BiGauss is the most precise and accurate diagnoses to measure 
even with relatively a poor sampling of the line profile. On the 
other hand, V^pan and W^pan need a good sampling of the line 
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This paper has been typeset from a Tj^/ ETeX file prepared by the 
author. 


Figure Cl. Top panel: Average time to measure the different spectroscopic 
diagnoses on a noiseless line profile as function of the sampling size (see 
text for details). Middle and bottom panels: Relative precision (middle 
panel) and accuracy (bottom panel) of the asymmetry diagnoses on a noise¬ 
less blended line profile as function of the sampling size. This assumes a 
line profile with a FWHM of lOkms^*. The RV and FWHM are not dis¬ 
played here since they are at the level of the numerical noise. 


profile to be measured precisely and accurately. This behaviour 
could be explained by the fact that the BiGauss diagnosis consists 
in fitting a continuous function, while Vspan and W^pa^ consist in 
fitting two discontinuous functions. BIS and Vasy follow the same 
precision and accuracy trend, and are dominated by the numerical 
noise produced by the interpolation of the line profile. The sam¬ 
pling of the line profile needs to be increased by a factor of more 
than 10 for BIS and Vasy to reach the same precision and accuracy 
as BiGauss. Given the computation time, the relative precision and 
accuracy, the sensitivity of the diagnoses to blends (Fig. 1) and to 
noise (see Section. 2.7), BiGauss is the most efficient spectroscopic 
diagnosis to constrain blended contamination. 
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Table Dl. Free parameters of the HD16702 analysis with their prior and posterior distributions. For the posterior we report only the median value and its 
68.3% uncertainty. ^{a,b) represents a Uniform prior between a and h; N(jd,cr^) represents a Normal distribution with a mean of fi and a width of cr^; 
Ti, b) represents a Normal distribution with a mean of /i and a width of cr^ and limited by a Uniform distribution between a and b. 


Parameter 

Prior 

Posterior 

Primary star 

Eifective temperature Tgff [K] 

^5908; 25) 

5909 ± 26 

Surface gravity log g [g.cm“^] 

yV(4.46;0.03) 

4.45 ± 0.03 

Iron abundance [Fe/H] [dex] 

M-0.12;0.03) 

-0.12 ±0.03 

Equatorial velocity 17 sin i [kms“^] 

AAm(2.56;1;0;20) 

3.42 ± 0.04 

Systemic radial velocity y [kms“^] 

<W(0; 10) 

4.576 ± 0.004 

Secondary star 

Initial mass M2 [M©] 

A^^/(0.55;0.14;0.1;20) 

0.35 ± 0.03 

Equatorial velocity r/sini* 2 [kms“^] 

15) 

6.4 ±2.1 

Binary orbit 

Orbital period P [d] 

^72.8322; 0.2) 

72.839 ± 0.003 

Periastron epoch Tp [BID - 2450000] 

A'(4983.45; 12) 

4983.20 ±0.14 

Eccentricity e 

ttiO-, 1) 

0.138 ±0.002 

Argument of periastron co [°] 

<W(0; 360) 

256.0 ±0.7 

Inclination i^rb [°] 

<7/(0; 90) 

9.4 ± 0.5 

Data-related parameters 

Radial velocity offset G2 - K5 [ m s ^ ] 

<W(-200;200) 

-1 ±5 

Radial velocity jitter G2 [ m s“ * ] 

-7/(0; 1000) 

15 ±3 

Radial velocity jitter K5 [ m s“ ^ ] 

-7/(0; 1000) 

15 ±3 

BiGauss offset G2 JbCgi 

‘77(-200;200) 

-1 ±3 

BiGauss offset K5 [ ni s“^] 

‘77(-200;200) 

-71 ±3 

BiGauss jitter G2 [ms“^] 

<77(0; 1000) 

3±3 

BiGauss jitter K5 ^ ] 

<77(0; 1000) 

4±4 

FWHM offset G2 - K5 Apw [ms'M 

<77(-200;200) 

7±6 

FWHM jitter G2 cr [ m s“ ‘ ] 

<77(0; 1000) 

8±6 

FWHM jitter K5 cr[ m s“ ^ 

<77(0; 1000) 

9±6 



